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Abstract. First part of the paper is devoted to applying the confinement mechanism proposed earlier by 
the author to estimate the possible parameters of the confining SU(3)-gluonic field in neutral kaons. The 
estimates obtained are consistent with the widths of the electromagnetic decays K° , K° — > 27 too. The 
corresponding estimates of the gluon concentrations, electric and magnetic colour field strengths are also 
adduced for the mentioned field at the scales of the mesons under consideration. The second part of the 
paper takes into account the results obtained previously by the author to estimate the purely gluonic 
contribution to the masses of all the mesons of pseudoscalar nonet and also to consider a possible relation 
with a phenomenological string-like picture of confinement. Finally, the problem of masses in particle 
physics is shortly discussed within the framework of approach to the chiral symmetry breaking in quantum 
chromodynamics (QCD) proposed recently by the author. 

PACS. 12.38.-t Quantum chromodynamics - 12. 38. Aw General properties of QCD (dynamics, confine- 
ment, etc.) - 14.40.Aq Pi, K, and eta mesons 
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1 Introduction 

The present paper to some degree summarizes our previ- 
ous ones on studying the pseudoscalar meson nonet within 
the framework of the (quark) confinement mechanism pro- 
posed earlier by the author. Global strategy may consist in 
reconsidering the whole spectroscopy of both mesons and 
baryons from the positions of the mentioned mechanism 
so exploring the pseudoscalar meson nonet is in essence 
only the first step in the given direction. But such a study 
involves concretizing the confinement mechanism itself in 
the generally accepted physical terms and when applying 
the mechanism to concrete hadrons there always arise new 
physical possibilities of interpreting the results obtained 
which further enriches the mechanism from physical point 
of view. In other words, the proposed confinement mecha- 
nism should be continuously modified and improved from 
physical positions and the best way of doing so is to study 
concrete hadrons with its help. Let us now shortly outline 
the main features of the approach suggested. 

In [UElEj f° r the Dirac- Yang-Mills system derived from 
QCD-Lagrangian an unique family of compatible nonper- 
turbative solutions was found and explored, which could 
pretend to decsribing confinement of two quarks. The ap- 
plications of the family to the description of both the 
heavy quarkonia spectra [H|5] and a number of properties 
of pions, kaons, 77- and ^'-mesons [5, 6,7,8,10] showed that 
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the confinement mechanism is qualitatively the same for 
both light mesons and heavy quarkonia. At this moment 
it can be decribed in the following way. 

The next main physical reasons underlie linear confine- 
ment in the mechanism under discussion. The first one is 
that gluon exchange between quarks is realized with the 
propagator different from the photon-like one, and exis- 
tence and form of such a propagator is a direct consequence 
of the unique confining nonperturbative solutions of the 
Yang-Mills equations [2J|3]. The second reason is that, ow- 
ing to the structure of the mentioned propagator, quarks 
mainly emit and interchange the soft gluons so the gluon 
condensate (a classical gluon field) between quarks basi- 
cally consists of soft gluons (for more details see Refs. [5J 
[3]) but, because of the fact that any gluon also emits glu- 
ons (still softer), the corresponding gluon concentrations 
rapidly become huge and form a linear confining magnetic 
colour field of enormous strengths, which leads to confine- 
ment of quarks. This is by virtue of the fact that just the 
magnetic part of the mentioned propagator is responsi- 
ble for a larger portion of gluon concentrations at large 
distances since the magnetic part has stronger infrared 
singularities than the electric one. In the circumstances 
physically the nonlinearity of the Yang-Mills equations ef- 
fectively vanishes so the latter possess the unique non- 
perturbative confining solutions of the Abelian-like form 
(with the values in Cartan subalgebra of SU(3)-Lie alge- 
bra) [21[3] which describe the gluon condensate under con- 
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sideration. Moreover, since the overwhelming majority of 
gluons is soft they cannot leave the hadron (meson) until 
some gluons obtain additional energy (due to an external 
reason) to rush out. So we also deal with the confinement 
of gluons. 

The approach under discussion equips us with the ex- 
plicit wave functions for every two quarks (meson or quarko- 
nium). The wave functions are parametrized by a set of 
real constants aj , bj , Bj describing the mentioned nonper- 
turbative confining SU(3)-gluonic field (the gluon conden- 
sate) and they are nonperturbative modulo square inte- 
grable solutions of the Dirac equation in the above confin- 
ing SU(3)-field and also depend on fio, the reduced mass 
of the current masses of quarks forming meson. It is clear 
that under the given approach just constants a 3 , bj , Bj , /xq 
determine all properties of any meson (quarkonium) , i. e., 
the approach directly appeals to quark and gluonic degrees 
of freedom as should be according to the first principles of 
QCD. Also it is clear that the constants mentioned should 
be extracted from experimental data. 

Such a program has been to a certain extent advanced 
in Rcfs. [4,5,6,7,8,9, 10] . Under the circumstances one aim 
of the present paper is to complete obtaining estimates 
for a j , bj , Bj for the nonet of pseudoscalar mesons and we 
shall here consider neutral kaons K°, K°. Another aim is 
to a certain degree to analyse some physical conclusions 
that can be obtained by considering the whole nonet from 
positions of the approach suggested. 

Of course, when conducting our considerations we shall 
rely on the standard quark model (SQM) based on SU(3)- 
flavor symmetry (see, e. g., so in accordance with 

SQM K°,K° = ds, ds respectively. 

Section 2 contains a survey of main relations underly- 
ing description of any mesons (quarkonia) in our approach. 
Section 3 gives estimates for parameters of the confining 
SU(3)-gluonic field for neutral kaons and also contains a 
discussion about whether the obtained estimates might 
also be consistent with the widths of 2-photon decays K°, 
K° — > 2j. Section 4 employs the obtained parameters of 
SU(3)-gluonic field to get the corresponding estimates for 
such characteristics of the mentioned field as gluon con- 
centrations, electric and magnetic colour field strengths 
at the scales of the mesons in question while Section 5 
deals with discussion about chiral limit for the nonet of 
pseudoscalar mesons. Section 6 takes into account the re- 
sults obtained previously by the author to consider a pos- 
sible relation with a phenomenological string-like picture 
of confinement. In section 7 the problem of masses in par- 
ticle physics is shortly discussed within the framework of 
approach to the chiral symmetry breaking in QCD pro- 
posed recently by the author. Section 8 is devoted to the 
concluding remarks. 

Appendices A and B contain the detailed description 
of main building blocks for meson wave functions in the 
approach under discussion, respectively: eigenspinors of 
the Euclidean Dirac operator on two-sphere 8 2 and radial 
parts for the modulo square integrable solutions of Dirac 
equation in the confining SU(3)-Yang-Mills field. At last, 
Appendix C supplements Section 2 with a proof of the 



uniqueness theorem from that Section in the case of SU(3)- 
Yang-Mills equations. 

Further we shall deal with the metric of the flat Minkowski 
spacetime M that we write down (using the ordinary set 
of local spherical coordinates r, ip for the spatial part) 
in the form 

ds 2 = g^dx^ l ®dx v = dt 2 -dr 2 -r 2 (d$ 2 + sin 2 tfdip 2 ), (1) 

so we have |<5| = | det(g A1 , y )| — (r 2 sin-#) 2 and < r < oo, 

< 1? < 7T, < <p < 2?r. 

Throughout the paper we employ the Heaviside-Lorentz 
system of units with h = c = 1, unless explicitly stated 
otherwise, so the gauge coupling constant g and the strong 
coupling constant a s are connected by the relation g 2 /(Air) = 
a s . 

When calculating we apply the relations 1 GeV -1 « 
0.1973269679 fm , 1 s" 1 w 0.658211915 x 10~ 24 GeV , 

1 V/m « 0.2309956375 xlO" 23 GeV 2 , 1 T = 47rxlO- 7 H/mx 

1 A/m sa 0.6925075988 x 10" 15 GeV 2 . 

Finally, for the necessary estimates we shall employ the 
Too-component (volumetric energy density ) of the energy- 
momentum tensor for a SU(3)-Yang-Mills field which should 
be written in the chosen system of units in the form 

V - -F; a F? p g a ? + i^ 5 Q VV • ( 2 ) 

2 Survey of main relations 

2.1 The confining SU(3)-gluonic field and meson wave 
functions 

As was mentioned above, our considerations shall be based 
on the unique family of compatible nonperturbative solu- 
tions for the Dirac- Yang-Mills system (derived from QCD- 
Lagrangian) studied at the whole length in Refs. [TJ[SJ[3]. 
Referring for more details to those references, let us briefly 
decribe and specify only the relations necessary to us in 
the present paper. 

One part of the mentioned family is presented by the 
unique nonperturbative confining solution of the SU(3)- 
Yang-Mills equations for the gluonic field A = A^dx 11 = 
A^Xadx^ (A a are the known Gell-Mann matrices, /i = 
t, r, #, (p, a = 1, 8) and looks as follows 

A xt = Al+±=A* = -^+A^A 2t ee ~^+^ = -y+^b, 

A 3t ^-^A! = ^-(A 1 + A 2 ), 
A lv = Al + -^=A% = b 1 r+B ll A 2v = -A s 9 +—Al = b 2 r+B 2 , 

A 3v = -^=A 8 v = -(b 1+ b 2 )r-(B 1+ B 2 ) (3) 

with the real constants aj , Aj , bj , Bj parametrizing the 
family. 
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The word unique should be understood in the strict 
mathematical sense. In fact in Ref. [2] the following theo- 
rem was proved (see also Appendix C): 

The unique exact spherically symmetric (nonperturba- 
tive) confining solutions (depending only on r and r _1 ) of 
SU (3) -Yang-Mills equations in Minkowski spacetime con- 
sist of the family of (3). 

It should be noted that solution (3) was found early 
in Ref. [I] but its uniqueness was proved just in Ref. [2] 
(see also Ref. [3])- Besides, in Ref. [5] (see also Ref. [5]) 
it was shown that the above unique confining solutions 
(3) satisfy the so-called Wilson confinement criterion [11] . 
Up to now nobody contested the above results so if we 
want to describe interaction between quarks by spheri- 
cally symmetric SU(3)-fields then they can be only those 
from the above theorem. On the other hand, the desirabil- 
ity of spherically symmetric (colour) interaction between 
quarks at all distances naturally follows from analysing 
the pp-collisions (see, e.g., Ref. [T3]) where one observes 
a Coulomb-like potential in events which can be identi- 
fied with scattering quarks on each other, i.e., actually at 
small distances one observes the Coulomb-like part of so- 
lution (3). Under this situation, a natural assumption will 
be that the quark interaction remains spherically symmet- 
ric at large distances too but then, if trying to extend the 
Coulomb-like part to large distances in a spherically sym- 
metric way, we shall inevitably come to the solution (3) 
in virtue of the above theorem. 

Now one should say that the similar unique confining 
solutions exist for all semisimple and non-semisimple com- 
pact Lie groups, in particular, for SU(iV) with N > 2 and 
U(iV) with N > 1 [213]. Explicit form of solutions, e.g., 
for SU(iV) with N = 2,4 can be found in Ref. [5] but it 
should be emphasized that components linear in r always 
represent the magnetic (colour) field in all the mentioned 
solutions. Especially, the case of the U(l)-group is inter- 
esting which corresponds to usual electrodynamics. Under 
this situation, as was pointed out in Refs. [2]|3], there is an 
interesting possibility of indirect experimental verification 
of the confinement mechanism under discussion. Indeed 
the confining solutions of Maxwell equations for classical 
electrodynamics point out the confinement phase could be 
in electrodynamics as well. Though there exist no elemen- 
tary charged particles generating a constant magnetic field 
linear in r, the distance from particle, after all, if it could 
generate this elecromagnetic field configuration in labora- 
tory then one might study motion of the charged particles 
in that field. The confining properties of the mentioned 
field should be displayed at classical level too but the exact 
behaviour of particles in this field requires certain analysis 
of the corresponding classical equations of motion. Such 
a program has been recently realized in Ref. [29] . Motion 
of a charged (classical) particle was studied in the field 
representing magnetic part of the mentioned solution of 
Maxwell equations and it was shown that one deals with 
the full classical confinement of the charged particle in 
such a field: under any initial conditions the particle mo- 
tion is accomplished within a finite region of space so that 
the particle trajectory is near magnetic field lines while 



the latter are compact manifolds (circles). Those results 
might be useful in thermonuclear plasma physics (for more 
details see [22]). 

As has been repeatedly explained in Refs. [2j[3l|4l[5], 
parameters Ai t % of solution (3) are inessential for physics 
in question and we can consider A\ = A 2 = 0. Also, as has 
been repeatedly discussed by us earlier (see, e. g., Refs. 
[2Mj), from the above form it is clear that the solution 
(3) is a configuration describing the electric Coulomb-like 
colour field (components A t ' ) and the magnetic colour 
field linear in r (components A^ 8 ) and we wrote down 
the solution (3) in the combinations that are just needed 
to insert into the corresponding Dirac equation. 

Another part of the family represents the meson wave 
functions and is given by the unique nonperturbative mod- 
ulo square integrable solutions of the mentioned Dirac 
equation in the confining SU(3)-field of (3) ^ = (tfq, ^3) 
with the four-dimensional Dirac spinors tyj representing 
the jth colour component of the meson, so W may de- 
scribe the relative motion (relativistic bound states) of 
two quarks in mesons and is at j — 1,2,3 (with Pauli 
matrix o~i) 

3 3 \ i^2(r)oi*i(0,¥O j w 

with the 2D eigenspinor <Pj = (^, 3l ^j °f t ne Euclidean 

Dirac operator T>o on the unit sphere S 2 , while the coor- 
dinate r stands for the distance between quarks. 

In this situation, if a meson is composed of quarks q\^ 
with different flavours then the energy spectrum of the 
meson will be given by e = m qi + m q2 +co with the current 
quark masses m qk ( rest energies) of the corresponding 
quarks and an interaction energy oj. On the other hand at 
3 = 1,2,3 

uij = uijirijjj^j) = 

A j g 2 a j b j ± (n 3 - + oij)J (n| + Znjccj + + g 2 b^(n^ + 2n$ocj) 

nj + 2n jaj + A] 

(5) 

with the gauge coupling constant g while /io is a mass pa- 
rameter and one should consider it to be the reduced mass 
which is equal to m qi m q2 / (m qi + m q2 ) with the current 
quark masses m qk ( rest energies) of the corresponding 
quarks forming a meson (quarkonium) , 03 = — [a\ + 02), 
b 3 = - (pi + 62), B 3 = -{B x + B 2 ), Aj = Xj - gB 0l 

ctj = yjAj - g 2 a 2 , rij = 0, 1,2, while Aj = ±(lj + 1) 

are the eigenvalues of Euclidean Dirac operator T>q on a 
unit sphere with lj = 0, 1, 2, .... 

In line with the above we should have lo = u>\ = L02 = 
u>3 in energy spectrum e = m qi + m q2 + u for any meson 
(quarkonium) and this at once imposes two conditions on 
parameters aj , bj , Bj when choosing some experimental 
value for e at the given current quark masses m qi , m q2 . 

The general form of the radial parts of (4) can be 
found, e.g., in Appendix B and within the given paper 
we need only the radial parts of (4) at nj = (the ground 
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state) that are 

Fji = cj'jr" . (i !J j-)- Pj = gbj + ft, 

F j2 = 'C,Q ,r- e~^ r (l + ^-j, Qj = Mo - Wj (6) 

with j3j = ^J/Iq — uj 2 + g 2 b 2 , while Cj is determined from 

the normalization condition f£°(\Fji\ 2 + \Fj2\ 2 )dr = |. 
The corresponding eigenspinors of (4) with A = ±1 (/ = 0) 
are 

A = -1 : * = % ( 6 %) e^ 2 ,or<Z> = £ ( e % ) e~W 

\ = l:$=^( e ~-*) e^ 2 , or* = § ( "C* ) e^/ 2 
2 V e 4 - / 2 V e ~ / 

(7) 

with the coefficient C = I/V^tt (for more details, see Ap- 
pendix A). 

2.2 Singularities of solutions 

As is seen from (3), the solutions in question have singular- 
ities: electric part contains the Coulomb-like singularities 
while we can rewrite the magnetic part in terms of differ- 
ential 1-forms as (bjr + Bj)dip, j =1, 2, 3, and then pass 
on to Cartesian coordinates employing the relations 

ox oy 

which entails 

n t~> \ i (hr + BAy , (bir + BAx , 

{b,r + Bj)d<p = — ^ ^-dx + 3 ^—dy , 

v J 3> f x 2 + y 2 x 2 + y 2 

wherefrom it is obvious that the colour magnetic field of 
(3) has the singularities on the z-axis. 

It should be noted that an analysis of singularites of 
YM-potentials requires both mathematical and physical 
considerations and, in general, is different from classical 
and quantum-mechanical point of view. In our recent pa- 
per [29] we gave some analysis of those singularities from 
classical point of view. 

But let us now note the following. In classical and 
quantum electrodynamics (QED) it is well known (see e. 
g. [16]) that the notion of classical electromagnetic field (a 
photon condensate) generated by a charged particle is ap- 
plicable only at distances much greater than the Compton 
wavelength A c = l/m for the given particle with mass m. 
Within the QCD framework the parameter Aqcd plays a 
similar part (see, e.g., [IS]) - Namely, the notion of classical 
SU(3)-gluonic field ( a gluon condensate) is not applicable 
at the distances much less than 1/ 'Aqcd- 

In this situation, the known singularity of the Coulomb 
potential in QED # = a/r at r = makes the purely 
mathematical sense since from the point of view of QED 



the photon condensate [huge number of (vitrual) photons] 
described by <P exists only at r >> \ c while at r < X c one 
may only speak about single photons rather than about 
condensate, i.e., the field in classical sense. The same holds 
true, e.g., for magnetic field of a uniformly moving charge 
where its strength H ~ v x r/r 3 , v is charge velocity. 

The colour magnetic field (3) under consideration has 
also the singularities on the z-axis so its formal math- 
ematical definition domain is the manifold R 3 \{z} with 
the z-axis discarded rather than the manifold R 3 . But we 
should not forget that we do not need only the appro- 
priate solutions of YM-equations to describe confinement. 
The YM-equations are only a part of the Dirac-YM sys- 
tem derived from QCD-lagrangian by standard prescrip- 
tion. The wave functions of hadrons (at any rate, mesons) 
are given by the modulo square integrable solutions of 
the Dirac equation (which is the second important part 
of the above Dirac-YM system) in the field (3). But the 
jth colour component of wave function of two quarks (me- 
son) (see (6)) in such a field behaves as ipj ~ r aj e~ 9 ^ b ^ r , 
[otj > 0), at \bj\ — > oo with bj characterising the linear 
colour magnetic field of solution (3), r is distance between 
quarks, j =1, 2, 3, 63 = —{b\ + 62)- he., typical size of 
hadron is r ~ l/(g\bj\) — > and we deal just with con- 
finement and besides we can see that those wave functions 
have no singularities along the z-axis and are well defined 
there, i.e. the wave functions are well defined already on 
the whole R 3 . But just the wave functions are needed to 
calculate miscellaneous characteristics of mesons (masses, 
radii and so on). So it is clear that at such computations 
the singularities of YM-potentials in questions have no in- 
fluence on physical results, i.e., those singularities are not 
physical ones. 

Physically this may mean that at large distances non- 
linearity of the Yang-Mills equations effectively vanishes 
so the latter possess the unique spherically symmetric non- 
perturbative confining solutions (3) (formally defined on 
R 3 \{z}) of the Abelian-like form (with the values in Car- 
tan subalgebra of SU(3)-Lie algebra) which describe the 
gluon condensate ( a classical gluon field) leading to the 
confinement. 

The situation is practically the same as for the hydro- 
gen atom or positronium: the wave functions of those sys- 
tems (see any textbook on quantum mechanics) are well 
defined at r — so the known singularity in the Coulomb 
potential at r — is also unphysical one, as said above. So 
modelling those singularities by some ^-functions (which 
is possible, as can easy show) makes no physical sense from 
the quantum-mechanical point of view and it can give just 
a suitable method for exploring some classical problems 
which is done, e.g., in many courses of classical electro- 
dynamics while the notion of classical chromodynamics in 
fact makes no sense: we can never generate the classical 
SU(3)-Yang-Mills field and classical coloured charged par- 
ticles at macroscopic scales. 

On the other hand, just quantum considerations can 
lead to one more point of view on the problem of singu- 
larity along z-axis of magnetic part for solution (3) and it 
is presented in Section 6 of the paper. 
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To summarize, solutions (3) are the unique spherically- 
symmetric solutions of YM-equations, though mathemat- 
ically being formally defined on M 3 \{z}, but with the un- 
physical singularities on the z-axis which are inessential 
from quantum-mechanical point of view. 



2.3 Nonrelativistic and the weak coupling limits 

It is useful to specify the nonrelativistic limit (when c — > 
oo ) for the spectrum (5). For this one should replace g — > 
g/y/hc, a,j — > cij/\/~hc, bj — > bj^fhc, Bj — > Bjj\/Wc and, 
expanding (5) in z = 1/c, we shall get 



IT 



9 2 a 2 



2V(n J + \\ 3 \y 



h(n 3 



//'/'/":/''',/■ 

fi3 (n . + | A .|)7 



z+0(z 2 ), (8) 



where f(rij, Xj) = 4X j n j (n* j +X})+^ (n$ T 6n|A? + A)). 

As is seen from (8), at c — > oo the contribution of linear 
magnetic colour field (parameters bj, Bj) to the spectrum 
really vanishes and the 

spectrum in essence becomes the purely nonrelativistic 
Coulomb one (modulo the rest energy). Also it is clear that 

when rij 



OO, UJj — > ± 



Mo 
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2 b 2 . At last, one should 



specify the weak coupling limit of (5), i.e., the case g ^ 0- 
As is not complicated to see from (5), ujj — > ±/^o when 

0. But then quantities f3j - 



Mo 



gib* 







g , ^^^^ — y f~u ~j 

and wave functions of (6) cease to be the modulo square 
integrable ones at g = 0, i.e., they cease to describe rel- 
ativistic bound states. Accordingly, this means that the 
equation (5) does not make physical meaning at g = 0. 

We may seemingly use (5) with various combinations 
of signes (±) before the second summand in numerators of 
(5) but, due to (8), it is reasonable to take all signs equal 
to plus which is our choice within the paper. Besides, as is 
not complicated to see, radial parts in the nonrelativistic 
limit have the behaviour of form Fji,Fj2 ~ r' J+1 , which 
allows one to call quantum number lj angular momentum 
for the jth colour component though angular momentum 
is not conserved in the field (3) 1,3 . So, for mesons under 
consideration we should put all lj = 0. 



2.4 Chiral limit 

There is one more interesting limit for relation (5) - the 
chiral one, i.e., the situation when m q i,m q 2 —> which 
entails /^o — > and (5) reduces to (at j = 1, 2, 3) 



Ajg 2 ajbj ± (n 3 + otj)g\bj 



j ) chiral 



2rijOij 



(9) 



which mathematically signifies that the Dirac equation 
in the field (3) possesses a nontrivial spectrum of bound 
states even for massless fermions. Physically this gives us 
a possible approach to the problem of chiral symmetry 
breaking in QCD [10] : in chirally symmetric world masses 
of mesons are fully determined by the confining SU(3)- 
gluonic field between (massless) quarks and not equal to 
zero. Accordingly chiral symmetry is a sufficiently rough 
approximation holding true only when neglecting the men- 
tioned SU(3)-gluonic field between quarks and no addi- 
tional mechanism of the spontaneous chiral symmetry break- 
ing connected to the so-called Goldstone bosons is re- 
quired. As a result, e.g., masses of mesons from pseu- 
doscalar nonet have a purely gluonic contribution and we 
shall consider it in section 5. 

One can note that for being the nonzero chiral limit 
of (5) the crucial role belongs to the colour magnetic field 
linear in r [parameters bi^ from solution (3)] inasmuch 
as chiral limit is equal exactly to zero when &1.2 = 0. On 
the contrary, when parameters ax,2 of the Coulomb colour 
electric part of solution (3) are equal to zero, the chiral 
limit may be nonzero at 61,2 7^ 0, as is seen from (9) except 
for the case rij = when both parts of SU(3)-gluonic field 
(3) are important for confinement and mass generation in 
chiral limit. 



2.5 Choice of quark masses and the gauge coupling 
constant 

Obviously, we should choose a few quantities that are the 
most important from the physical point of view to char- 
acterize mesons under consideration and then we should 
evaluate the given quantities within the framework of our 
approach. In the circumstances let us settle on the ground 
state energy (mass) of neutral kaons, the root-mean-square 
radius of them and the magnetic moment. All three mag- 
nitudes are essentially nonperturbative ones, and can be 
calculated only by nonperturbative techniques. 

Within the present paper we shall use relations (5) at 







so energy (mass) of mesons under considera- 



tion is given by /i = m<2 + m s +lu with u> — lu 3 (0, 0, Xj) for 
any j = 1,2,3 whereas 



g 2 a 1 b 1 «i^ g 2 a 2 b 2 , a 2 ^o 



Ax 

T, 



+ 



\A 3 



\i - m d 



\A 2 



(10) 



and, as a consequence, the corresponding meson wave func- 
tions of (4) are represented by (6) and (7). It is evident 
for employing the above relations we have to assign some 
values to quark masses and gauge coupling constant g. 
We take the current quark masses used in [MTIISIIWTU] 
and they are md — 5 MeV, m s — 107.5 MeV. Under the 
circumstances, the reduced mass /j,o of (5) will be equal 
to mdTn s / (m,d + m s ). As to the gauge coupling constant 
g — yJAitczs, it should be noted that recently some at- 
tempts have been made to generalize the standard for- 
mula for a s = a s (Q 2 ) = 12 7 r/[(33-2n / )ln(Q 2 //l 2 )] (n f is 
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number of quark flavours) holding true at the momentum 
transfer y/Q 2 — > oo to the whole interval < y/Q 2 < oo . 
If employing one such a generalization used in Refs. [H] 
which we have already discussed elsewhere (for more de- 
tails see [MTUMMin]) then (when fixing A = 0.234 GeV, 
ri/ = 3) we obtain g ps 5.290449085 necessary for our 
further computations at the mass scale of neutral kaons. 



2.6 Electric form factor and the root-mean-square 
radius 

The relations (4), (6) and (7) allow us to compute an 
electric formfactor of a meson as a function of the square 
of momentum transfer Q 2 in the form (for more details 



f(Q 2 



£ 



(2ft) 



2a, +1 



6a, 



iin [2aj arctan (y/\Q 2 \/(20j))] 

VWK^ - Q 2 Y 3 



(ii) 



which also entails the root-mean-square radius of the me- 
son (quarkonium) in the form 



<r>- 



Y.— 



3aj + 1 



6/3 



(12) 



that is in essence a radius of confinement. 
2.7 Magnetic moment 

Also it is not complicated to show with the help (4), (6) 
and (7) that the magnetic moments of mesons (quarkonia) 
with the wave functions of (4) (at lj = 0) are equal to zero 
[6.7,8,9,10 , as should be according to experimental data 

Though we can also evaluate the magnetic form factor 
F(Q 2 ) of meson (quarkonium) which is also a function of 
Q 2 (see Refs. [HIE] ) the latter will not be used in the given 
paper so we shall not dwell upon it. 

3 Estimates for parameters of SU(3)-gluonic 
field in neutral kaons 

3.1 Basic equations and numerical results 

Now we are able to estimate parameters a,j , bj , Bj of the 
confining SU(3)-field (3) for neutral kaons within frame- 
work of our approach. In this situation, we should consider 
(10) and (12) the system of equations which should be 
solved compatibly if taking fi = 497.648 MeV, = 5.0 
MeV, m s = 107.5 MeV and < r >« 0.560 fm in accor- 
dance with [13] . While computing for distinctness we take 
all eigenvalues Aj of the Euclidean Dirac operator T>q on 
the unit 2-sphere § 2 equal to 1 . The results of numerical 
compatible solving of equations (10) and (12) are adduced 
in Tables 1-2. 



3.2 Consistency with the widths of 2-photon decays 

K°, K° -> 2 7 

Let us consider whether the estimates of previous subsec- 
tion are consistent with the width of the electromagnetic 
2-photon decays K°, K° — > 27. Actually kinematic analy- 
sis based on Lorentz- and gauge invariances gives rise to 
the following expression for the width P of the electro- 
magnetic decay P — > 27 (where P stands for any meson 
from tt°, 77, 77', K°,K°, see, e.g., Ref. [TS]) 



r = \^al m g 2 Pll [i z 



(13) 



with the electromagnetic coupling constant a em =1/137.0359895 
and the P-meson mass \x while the information about 
strong interaction of quarks in P-meson is encoded in a 
decay constant <Jp 77 . Making replacement <?p 77 = fp/n 
we can reduce (13) to the form 



r 



™ 2 em vf'p 



(14) 



Now it should be noted that the only invariant which f P 
might depend on is Q 2 = [i 2 , i. e. we should find such a 
function F{Q 2 ) for that F(Q 2 = /j 2 ) = f P but F{Q 2 ) can- 
not be computed by perturbative techniques. It is obvious 
from the physical point of view that T(Q 2 ) should be con- 
nected with the electromagnetic properties of P-meson. 
As we have seen in Section 3, there are at least two suit- 
able functions for this aim - electric and magnetic form 
factors. But there exist no experimental consequences re- 
lated to a magnetic form factor at present whereas electric 
one to some extent determines, e. g., an effective size of 
meson (quarkonium) in the form < r > of (12). It is rea- 
sonable, therefore, to take T(Q 2 = fi 2 ) = Af(Q 2 = fi 2 ) 
with some constant A and the electric form factor / of 
(11) for the sought relation. In this situation, we obtain 
an additional equation imposed on parameters of the con- 
fining SU(3)-gluonic field in P-meson which has been used 
in Refs. [HlE] to estimate the mentioned parameters in tt°- 
and 77-mesons. As a result, using (11) in the case of neutral 
kaons, we come from (14) to relation 



r 



A 



3 



1 sin(2a J arctana; :) ) 



' 6a, 



(1 



0.209 x 10- 10 eV, K° s 
0.696 x 10- n eV, K° L 



j > 

mode, 
mode 



(15) 



with xj = /V(2/3j), fJ- = 497.648 MeV and we used widths 
r 7 w_0.209 x 10- 10 eV, A 7 « 0.696 x 10" 11 eV for decays 
K°,K° — > 27, respectively, for Kg - and if^-modes fol- 
lowing the notation from Ref. |13j . In the circumstances, 
we can employ the results of Table 1 and compute the left- 
hand side of (15) which entails the corresponding values 
A « 0.3263 x 10~ 7 and A w 0.1884 x 10~ 7 . Consequently, 
we draw the conclusion that parameters of the confining 
SU(3)-gluonic field in neutral kaons from Table 1 might be 
consistent with P7 and A7 while smallness of constants A 
indicates the electromagnetic properties of neutral kaons 
to be inessential. 
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Table 1. Gauge coupling constant, reduced mass /Jo and parameters of the confining SU(3)-gluonic field for neutral kaons 



Particle 


9 


Mo (MeV) 


di 


«2 


bi (GeV) 


b 2 (GeV) 


Si 


B 2 


K°,K° - ds, ds 


5.29045 


4.77778 


0.102484 


-0.198658 


0.385250 


-0.130208 


-0.360 


-0.170 



Table 2. Theoretical and experimental mass and radius of neutral kaons 



Particle 


Theorct. fi (McV) 


Experim. it (MeV) 


Theoret. < r > (fm) 


Experim. < r > (fin) 


K°,K° - ds., ds 


fi - m d + rn a + uij(0, 0, 1) = 497.648 


497.648 


0.550510 


0.560 



4 Estimates of gluon concentrations, electric 
and magnetic colour field strengths 

Now let us recall that, according to Refs. 3.5 . one can 
confront the field (3) with the Too-component (the volu- 
metric energy density of the SU(3)-gluonic field) of the 
energy-momentum tensor (2) so that 



T, 



oo 



T a = 



E z + H 2 



aia 2 



b\ + bib 2 



r 2 sin 2 



A 

,,4 



B 



r 2 sin $ 
(16) 

with electric E and magnetic H colour field strengths and 
with real A > 0, B > 0. One can also introduce magnetic 
colour induction B — {An x lCP 7 H/m) H, where H in 
A/m. 

To estimate the gluon concentrations we can employ 
(16) and, taking the quantity uj = T, the full decay width 
of a meson, for the characteristic frequency of gluons we 
obtain the sought characteristic concentration n in the 
form 

Too 



r 



so we can rewrite (16) in the form Too = T§q u1 + T uo n con- 
forming to the contributions from the Coulomb and linear 
parts of the solution (3). This entails the corresponding 
split of n from (17) as n = n cou i + ny ln . 

The parameters of Table 1 were employed when com- 
puting and for simplicity we put sini? = 1 in (16). There 
was also used the following present-day full decay widths 
of mesons under consideration [13 : r = 1/r with the life 
times t = 0.8953 x 1(T 10 s (tfg-mode) , 5.18 x 1CT 8 s 
(i"T^-mode), respectively, whereas the Bohr radius clq = 
0.529177249 • 10 5 fm [T3]. 

Table 3 contains the numerical results for n cou i, nu n , 
7i, E, H, B for the mesons under discussion. 



4.1 Concluding remarks 

As is seen from Table 3, at the characteristic scales of 
neutral kaons the gluon concentrations are huge and the 
corresponding fields (electric and magnetic colour ones) 
can be considered to be the classical ones with enormous 



strenghts. The part n cou i of gluon concentration n con- 
nected with the Coulomb electric colour field is decreas- 
ing faster than nn n , the part of n related to the linear 
magnetic colour field, and at large distances nu n becomes 
dominant. It should be emphasized that in fact the gluon 
concentrations are much greater than the estimates given 
in Table 3 because the latter are the estimates for max- 
imal possible gluon frequencies, i.e. for maximal possible 
gluon impulses (under the concrete situation of neutral 
kaons). As was mentioned in section 1, the overwhelm- 
ing majority of gluons between quarks is soft, i. e., with 
frequencies much less than r = 1/r with the life time 
r for Kg or Kl, so the corresponding concentrations are 
much greater than those in table 3. The given picture is in 
concordance with the one obtained in 4,5,6,7,9,10 . As 
a result, the confinement mechanism developed in [TJ[2l[3] 
and described early in section 1 is also confirmed by the 
considerations of the present paper. 

It should be noted, however, that our results are of a 
preliminary character which is readily apparent, for ex- 
ample, from the fact that the current quark masses (as 
well as the gauge coupling constant g) used in computa- 
tion are known only within the certain limits, and we can 
expect similar limits for the magnitudes discussed in the 
paper so it is neccesary for further specification of the pa- 
rameters for the confining SU(3)-gluonic field in neutral 
kaons which can be obtained, for instance, by calculating 
the widths of decays Kg —> 7r + 7r~ or K^ 3tt° with the 
help of wave functions discussed above and in [fjlfTU], We 
hope to continue analysing the given problems elsewhere. 



5 Chiral limit for pseudoscalar nonet 

Having obtained estimates for neutral kaons in previous 
sections we can state that at the given moment we have 
such estimates for all the members of pseudoscalar nonet 
if taking into account the results of [6l[7ll8lfT0] . Under the 
circumstances we can return to the chiral symmetry break- 
ing problem in QCD whose possible resolution within the 
framework of the above confinement mechanism has been 
discussed in [10 . As was mentioned in section 2, merits 
of case consists in that the Dirac equation in the field (3) 
possesses a nontrivial spectrum of bound states even for 
massless fermions [(see relation (9)]. As a result, mass of 
any meson remains nonzero in chiral limit when masses 
of quarks m q — > and meson masses will only be ex- 
pressed through the parameters of the confining SU(3)- 
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Table 3. Gluon concentrations, electric and magnetic colour field strengths in neutral kaons 



Kg-modc: r =< r >= 0.550510 fm 



r (fm) 


"ooul fm" 3 ) 


n 


lin (™ - 


3 




n 


(„-3) 






E (V/m) 






H (A/m) 




B 


(T) 






O.lro 


0.285346 X 10 65 





336488 


X 


10 63 





288711 




10 65 


0.957080 


x 


10 24 


0.139807 X 


10 22 





175687 


x 


10 16 


r 


0.285346 X 10 61 





336488 


X 


10 61 





621834 




10 61 


0.957080 


X 


10 22 


0.139807 X 


10 21 





175687 


X 


10 15 


1.0 


0.262079 X 10 60 





101976 




10 61 





128184 




10 61 


0.290054 


X 


10 22 


0.769654 X 


10 20 





967175 




10 14 


10r 


0.285346 X 10 57 





336488 




10 59 





339342 




10 59 


0.957080 




10 20 


0.139807 X 


10 20 





175687 




10 14 


"0 


0.334217 x 10 41 





364165 




10 51 





364165 




10 51 


0.103580 




10 13 


0.145443 X 


10 16 





182770 




10 10 


K£-mode: 


r =< r >= 0.550510 fm 




































r (fm) 


-com C-- 3 ) 




lin Crn" 


3 






(m-3) 






E (V/m) 






H (A/m) 




B 


(T) 






0.1r 


0.165095 x 10 6S 





194684 


X 


10 66 





167042 


X 


10 68 


0.957080 


X 


10 24 


0.139807 X 


10 22 





175687 


X 


10 16 


r 


0.165095 x 10 64 





194684 


X 


10 64 





359779 


X 


10 64 


0.957080 


X 


10 22 


0.139807 X 


10 21 





175687 




10 15 


1.0 


0.151633 x 10 63 


(") 


590013 


X 


10 63 





741646 


X 


10 63 


0.290054 


X 


10 22 


0.769654 X 


10 20 





967175 


X 


10 14 


10r 


0.165095 x 10 60 





194684 


X 


10 62 





196335 


X 


10 62 


0.957080 


X 


10 20 


0.139807 X 


10 20 





175687 


X 


10 14 


a 


0.193370 X 10 44 





210697 


X 


10 54 





210697 


X 


10 54 


0.103580 




10 13 


0.145443 X 


10 16 





182770 


X 


10 10 



gluonic field of (3). This purely gluonic residual mass of 
meson should be interpreted as a gluonic contribution to 
the meson mass. 

Physically this gives us a possible approach to the 
problem of chiral symmetry breaking in QCD [10] : in chi- 
rally symmetric world masses of mesons are fully deter- 
mined by the confining SU(3)-gluonic field between (mass- 
less) quarks and are not equal to zero. Accordingly chiral 
symmetry is a sufficiently rough approximation holding 
true only when neglecting the mentioned SU(3)-gluonic 
field between quarks and no additional mechanism of the 
spontaneous chiral symmetry breaking connected to the 
so-called Goldstone bosons is required. Referring for more 
details to [TOJ, we can here only say that, e.g., masses 
of mesons from pseudoscalar nonet have a purely gluonic 
contribution and we may be interested in what part of the 
meson masses is obligatory to that contribution. There- 
fore, let us employ the results of both [6,7,8,10 and the 
present paper to estimate the mentioned contribution for 
all the members of pseudoscalar nonet. To pass on to ob- 
taining the sought estimates all the necessary parameters 
gained in [61I71IMTU] are gathered in table 4 where we took 
into account that in accordance with the standard quark 
model based on SU(3)-flavour symmetry (see, e.g., 13;) 
7T° = (uu — dd)/y/2 is a superposition of two quarkonia 
while r\ = (2ss — uu — dd) / v6 and rj = (uu + dd + ss)/ 
are the superpositions of three quarkonia so we have, re- 
spectively, two or three sets of parameters aj, bj, Bj for 
the corresponding particles. 

The current quark masses were taken with the same 
values as in |g)l7llg|IT0"] . i.e., m u = 2.25 MeV, m d = 5 MeV, 
m s = 107.5 MeV. 

Now we can note that according to our approach the 
mass of any meson from table 4 is given by relation [cf. 
(10)] 



m a 



m 



g 2 aih 
Ai - gB 1 



Mo 



v / (Ai-.gSi) 2 



g 2 a\ 



\Xx-gBx 



g 2 a 2 b2 
A2 - gB 2 



Mo 



v/(A 2 - gB 2 y 



g 2 a\ 



|A 2 - gB 2 



''a 



''12 



g 2 (ai +a 2 )(fri + b 2 ) 
A 3 + 5 (Bi + ^2) 



/'o 



y/[X 3 + ,g(Bi + B 2 )] 2 - g 2 {ai 
\h+g(B! + B 2 )\ 



Mo 



'91 '"92 



(18) 



and the eigenvalues of the Euclidean Dirac operator on a 
unit sphere Xj 
particles. 



T for 7T° and A, 



1 for the rest of 



In chiral limit m qi , m q2 



chiral 



we obtain 
g 2 ai&i g 2 a 2 b 2 



g 2 (ai 



Ai - gBi A 2 - gB 2 
+ a 2 ){b 1 +b 2 ) 



h + g{Bi + B 2 ) 



7^0. 



(19) 



We can see that in chiral limit the meson masses are com- 
pletely determined only by the parameters aj , bj , Bj of 
SU(3)-gluonic field between quarks, i.e. by interaction be- 
tween quarks, and those masses have the purely gluonic 
nature. So one can use the parameters g, cij, bj, Bj adduced 
in table 4 to compute (n) c hirai which in fact represents the 
sought gluonic contribution to the meson masses. The re- 
sults of computation are gathered in table 5. 

One can add to the results of table 5 that we could 
also calculate, e.g., the root-mean-square radii (12) of par- 
ticles under consideration in chiral limit which are well 
defined in this limit as well. For example, < r > c hirai~ 
0.673069 fm or 0.543223 fm, accordingly, for charged pions 
and kaons jTD] . I.e., those values only slightly differ from 
the present-day experimental values 0.672 fm or 0.560 fm 
|13j . The same remark also holds true for the decay con- 
stants for leptonic decays fp (P stands for charged pi- 
ons and kaons, see [TU] for more details). So, even in chi- 
rally symmetric world, e.g., the charged pions and kaons 
would have nonzero masses, the root-mean-square radii 
and decay constants fp for leptonic decays and all of those 
quantities would be determined only by SU(3)-gluonic in- 
teraction between massless quarks, i.e. they would have 
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Table 4. Gauge coupling constant, reduced mass /io and parameters of the confining SU(3)-gluonic field for pseudoscalar nonet. 



Particle 


9 


fl (MeV) 


ai 


a 2 


6i (GeV) 


62 (GeV) 


B 1 


B 2 


n - 

7T — UU 


6.10148 


1.125 


-0.0434737 


-0.00680835 


0.0848234 


0.0433136 


0.01 


-0.150 


7r n - dd 


6.10148 


2.50 


-0.0606679 


0.0251427 


0.0956303 


0.0648174 


0.1250 


-0.2450 


tt^ 1 - ud, ud 


6.09131 


1.55172 


0.0473002 


0.0118497 


0.178915 


-0.119290 


-0.230 


0.230 


K^ 1 — us, us 


5.30121 


2.20387 


0.167182 


-0.0557501 


0.120150 


0.131046 


-0.900 


0.290 


K°,K° - ds, ds 


5.29045 


4.77778 


0.102484 


-0.198658 


0.385250 


-0.130208 


-0.360 


-0.170 


7} — UU 


5.14836 


1.125 


-0.0328122 


0.179728 


0.194979 


0.119737 


0.255 


-0.010 


r/ — dd 


5.14836 


2.50 


0.147640 


-0.178707 


0.305728 


-0.119050 


-0.240 


-0.010 


7] — SS 


5.14836 


53.75 


-0.0141391 


-0.0806779 


0.252975 


-0.339250 


0.260 


-0.310 


7] - UU 


3.91476 


1.125 


0.218474 


-0.394718 


0.618419 


-0.280807 


-0.300 


-0.200 


7] - dd 


3.91476 


2.50 


0.351384 


-0.130858 


0.278983 


0.285548 


-0.150 


0.410 


7}' - SS 


3.91476 


53.75 


0.123645 


0.124633 


-0.226875 


0.588802 


0.410 


-0.160 



Table 5. Theoretical, experimental and chiral meson masses 



Particle 


Theoretical (MeV) 


Experimental (MeV) 


Chiral (MeV) 


Gluonic contribution (%) 


7T° — UU 


u = 2m u + cjj(0, 0, -1) = 134.976 


134.976 


129.495 


96.0 


7r° - dd 


u = 2m d + Wj(0, 0, -1) = 134.976 


134.976 


122.5 


90.8 


7r^ — ud. ud 


u = ro„ + m d + uij (0, 0, 1) = 139.570 


139.56995 


130.8 


93.7 


— us, us 


u = m u + m a + t»j(0, 0, 1) = 493.677 


493.677 


382.0 


77.4 


K°,K° - ds, ds 


u = m d + m a + Uj (0, 0, 1) = 497.648 


497.648 


380.569 


76.6 


rj — uu 


u = 2ra„ + w,(0, 0, 1) = 547.51 


547.51 


542.0 


99.0 


rj — dd 


u = 2m d + uj (0, 0, 1) = 547.51 


547.51 


535.4 


97.8 


rj — ss 


u = 2m s + to, (0. 0, 1) = 547.51 


547.51 


280.01 


51.1 


rj' — uu 


u = 2m u + UjiO, 0, 1) = 957.78 


957.78 


952.5 


99.4 


rj' — dd 


u = 2m d + uj(0, 0, 1) = 957.78 


957.78 


946.0 


98.8 


rj' — ss 


u = 2m 3 + ujj (0. 0. 1) = 957.78 


957.78 


691.5 


72.2 



a purely gluonic nature. Moreover, since gluons are ver- 
ily relativistic particles then the most part of masses for 
mesons of pseudoscalar nonet is conditioned by relativis- 
tic effects, as is seen from table 5. Further discussion of 
the proposed chiral symmetry breaking mechanism can be 
found in [TO] . 

6 A possible relation with a phenomenological 
string-like picture of quark confinement 

6.1 The confining potential and string tension 

The results adduced in section 5 allow us to shed some 
light on one more problem which has been touched upon in 
[5, 10]. As is known, for a long time up to now there exists 
the so-called string-like picture of quark confinement but 
only at qualitative phenomenological level (see, e. g., Ref. 
[12]). Up to now, however, it is unknown how such a pic- 
ture might be warranted from the point of view of QCD. 
Let us in short outline as our results for pseudoscalar 
nonet (based on and derived from QCD-Lagrangian di- 
rectly) naturally lead to possible justification of the men- 
tioned contruction. Thereto we note that one can calculate 
energy £ of gluon condensate conforming to solution (3) 
in a volume V through relation £ = f v Tqqt 2 sin ddrd'ddip 
with Tqo of (16) but one should take into account that 



classical Too has a singularity along z-axis (i9 = 0, tt) and 
we have to introduce some angle 1)0 so i)q < -d < tt — i}q. 
As well as in Ref. [5] , we may consider $0 to be a param- 
eter determining some cone d = $0 so the quark emits 
gluons outside of the cone. Now if there are two quarks 
Qi, Q2 and each of them emits gluons outside of its own 
cone i9 = $1,2 (see Figs. 1, 2) then we have soft gluons 
(as mentioned in section 1) in regions I, II and between 
quarks. 

Accordingly, we shall have some region V with gluon 
condensate between quarks Qi, Q2 and its vertical projec- 
tion is shown in Fig. 1. Another projection of V onto a 
plane perpendicular to the one of Fig. 1 is sketched out in 
Fig. 2. 

Then, as is clear from Fig. 1, for distance R between 
quarks we have R — R± sin^i + R2 sin ^2 and gluonic en- 
ergy between quarks will be equal to 

V{R)= { T 00 r 2 sin Mr d#d<p = 
Jv 

with constants A, B defined in (16). 



10 
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Fig. 1. Vertical projection of region with the gluon condensate energy between quarks. 



Fig. 2. Horizontal projection of region with the gluon condensate energy between quarks. 



To clarify a physical meaning of the quantities in 
Figs. 1, 2, let us recall an analogy with classical electrody- 
namics where is well known (see, e. g., |16] and Subsection 
2.2) that the notion of classical electromagnetic field (a 
photon condensate) generated by a charged particle is ap- 
plicable only at distances much greater than the Compton 
wavelength A c = 1/m for the given particle with mass m. 
Within the QCD framework the parameter Aqcd plays a 
similar part (see, e.g., Ref. [T3IIT7] ). Namely, the notion of 
classical SU(3)-gluonic field ( a gluon condensate) is not 
applicable at the distances much less than 1/ Aqcd- In ac- 
cordance with subsection 2.5 we took Aqcd — A — 0.234 
GeV which entails 1/A ~ 0.8433 fm so one may consider 
n,2 ~ 0.1< r > with the root-mean-square radius < r > 
for a meson. 

Under the circumstances, performing a simple integra- 
tion in (20) with employing the relations J <i0/sin0 = 
In tan. tan0/2 = sin0/(l + cos0) = (1 - cos0)/sin0, 
we shall without going into details (see also Ref. |2j) ob- 
tain 



,,/„ o n ^ tyiAcostfi A l + cos0,: 

V(Ri,R 2 ) = V -2^ □ h> 2(piBRi\sx — 

i=i 1 i=i 



where 



cos i9i 
(21) 



Vo = X>ch 



i=l 



E 



/ AtpiA cos §i 

\ n 



2ip. l Br l In 



1 + cos 0j 
1 — cos 0, 



For the sake of simplicity let us put R\ = R2, 0i = 
02 = 0o, <pi = fi = fo- Then R l = R 2 = i?/(2sin0 o ) 
and from (21) it follows 



V(.R) = Vo + - + kR 



(22) 



with a = -S^sin 20 o , k = 2^^ In ±±^|fa. 

We recognize the modeling confining potential in (22) 
which is often used when applying to meson and heavy 
quarkonia physics (see, e.g., [H]). We can, however, see 
that phenomenological parameters a, k, Vo of potential (22) 
are expressed through more fundamental parameters Oj, 
bj connected with the unique exact solution (3) of Yang- 
Mills equations describing confinement. One can notice 
that the quantity k (string tension) is usually related to 
the so-called Regge slope a' = l/(2nk) and in many if not 
all of the papers using potential approach it is accepted 
k w 0.18 GeV 2 (see, e. g., [IB]). 



6.2 Estimates of O , ipo for pseudoscalar nonet 

Under the situation we have the equation 



k = 2ipo~ 



B 



1 + COS 0o 

sin 0o 1 — cos 0q 



-In 



> 0.18 GeV 2 (23) 



with B = (b\ + bib 2 + b\)j2, so let us employ (23) to esti- 
mate 0o , <po if using the parameters adduced in table 4 for 
pseudoscalar nonet and also for the ground state of topo- 
nium rjt for that we use the parametrization from [9 with 
the values a x = 0.361253, a 2 = 0.339442, 61 = 48.9402 
GeV, b 2 = 76.7974 GeV for the parameters of solution 
(3). Results of computations are presented in table 6. 

If taking into account that only the values of 0o, fo 
between and 90° are of physical meaning and, accord- 
ing to Figs. 1, 2, the gluon configuration between quarks 
will be similar to a string-like one under the condition 
0o ^ > t/2, ipo — > 0, then we can see from table 6 that the 
characteristic transverse sizes D\ 2 of the gluon conden- 
sate between quarks in fact tend to zero only in the case 
of heavy quarks, i.e., only for heavy quarks the gluon con- 
figuration between them might practically transform into 
a string. As a result, there arises the string- like picture of 
quark confinement but the latter seems to be warranted 
enough only for heavy quarks. It should be emphasized 
that string tension k of (23) is determined just by param- 
eters &i, 2 of linear magnetic colour field from solution (3) 
which indirectly confirms the dominant role of the men- 
tioned field for confinement. 

We cannot, however, speak about potential V(i?) of 
(22) as describing some gluon configuration between quarks. 
It would be possible if the mentioned potential were a so- 
lution of Yang-Mills equations directly derived from QCD- 
Lagrangian since, from the QCD-point of view, any glu- 
onic field should be a solution of Yang-Mills equations (as 
well as any electromagnetic field is by definition always a 
solution of Maxwell equations) . 

In reality, as was shown in Refs. [2]|3] (see also Ap- 
pendix C), potential of form (22) cannot be a solution 
of the Yang-Mills equations if simultaneously o ^ 0,fc ^ 
0. Therefore, it is impossible to obtain compatible solu- 
tions of the Yang-Mills-Dirac (Pauli, Schrodinger) system 
when inserting potential of form (22) into Dirac (Pauli, 
Schrodinger) equation. So, we draw the conclusion (men- 
tioned as far back as in Refs. [1] and elaborated more in de- 
tail in Ref. [3]) that the potential approach seems to be in- 
consistent: it is not based on compatible nonperturbative 
solutions for the Dirac- Yang-Mills system derived from 
QCD-Lagrangian in contrast to our confinement mecha- 
nism. Actually potential approach for heavy quarkonia has 
been historically modeled on positronium theory. In the 
latter case, however, one uses the unique modulo square 
integrable solutions of Dirac (Schrodinger) equation in the 
Coulomb field [condensate of huge number of (virtual) 
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Table 6. Angular parameters determining the gluon condensate between quarks for pseudoscalar nonet and toponium ground 
state 





fin 






10° 


28 84° 




30° 


153 6° 


-TT — dd 


10° 


18 8° 




30° 


100 2° 


± A 
IX UCL > 


30° 


78 63° 




45° 


166 16° 


jV US ^ US 


45° 


87 36° 




60° 


171 68° 


K Q K° — ds da. 


60° 


70 6° 




70° 


118 0° 




45° 


54.7° 


— — i 


60° 


107 4° 


Id 


45° 


58 0° 


— - 1 


60° 


114.1° 


rj - ss 


60° 


87.2° 




70° 


145.8° 


rj — uu 


70° 


47.3° 




80° 


100.6° 


rj — dd 


70° 


56.9° 




80° 


121.1° 


rj — ss 


70° 


75.5° 




80° 


160.8° 


rjt - it 


60° 


(0.675 x 10" 3 )° 




80° 


(0.240 x 10 _2 )° 




88° 


(0.123 x 10 _1 )° 



photons], i. e., one employs the unique compatible non- 
perturbative solutions of the Maxwell-Dirac (Schrodinger) 
system directly derived from QED-Lagrangian to describe 
positronium (or hydrogen atom) spectrum. 

To summarize, from the point of view of our approach 
both potential and string- like pictures of confinement arise 
only as some effective models derived in a certain way 
from the more fundamental theory based on exact solution 
(3) of SU(3)-Yang-Mills equations. This conlusion is in 
concordance with the ones obtained in I51HU1. 



charge. If taking into account that the mass of Higgs bo- 
son tuh = s\v with a self-interaction constant A then it 
is clear that masses of all the abovementioned particles 
are proportional to to#, and, consequently, the discovery 
of Higgs boson will not completely resolve the puzzle of 
origin of masses in particle physics - the question will re- 
main where the mass ran comes from not speaking already 
about the nature of the above miscellaneous constants fi 
and A. 

At present, to our mind, one can single out two most 
promising approaches to a possible resolution for the men- 
tioned problems: technicolour theories and preon models. 
Under the circumstances let us shortly outline how both 
these directions might be estimated from the point of view 
of our confinement mechanism and the chiral symmetry 
breaking one based on the latter and discussed above and 
in [TO]. 

7.2 Technicolour theories 

Referring for more details concerning those models to both 
early references [TH] and modern status of them (see, e.g., 
ptlj l let us note the following. The main idea of acquir- 
ing masses, e.g., for and Z bosons, consists in that 
a new set of the so-called techniquarks is postulated at 
for the gauge bosons , Z where masses raw = ev/(2 sin^vit^e energy scale of order 1 TeV which interact with each 
mz = ev/(sm26w) with the so-called weak angle 9\y other through the technigluons and it makes the mass- 
so that sin 2 9w ~ 0.23 and e is the elementary electric less techipions exist as Goldstone bosons. The latter give 



7 Problem of masses in particle physics 
7.1 Preliminaries 

As is known |13) . the generally accepted standard model 
with one Higgs doublet asserts that the masses of funda- 
mental fermions (quarks and leptons) are acquired through 
the Higgs mechanism so for their masses m.j we obtain 
(without taking mixings into account) nij = /ju/v2, where 
the vacuum Higgs condensate v w 246 GeV and i stands 
for quark and lepton flavours. But little is known about 
the coupling constants /, and much may be elucidated 
only with discovering Higgs bosons. The same holds true 
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masses to W and Z after spontaneous symmetry break- 
ing. It should be noted, however, those massless technip- 
ions appear as a result of violating chiral symmetry con- 
nected with technicolour QCD on the analogy with chiral 
symmetry breaking in usual QCD. But, as we have dis- 
cussed in [TU] and in section 5, the hypothetical mecha- 
nism for chiral symmetry breaking in standard QCD with 
appearance of Goldstone bosons (pions) seems to fail be- 
cause of pions can never be massless inasmuch as they have 
nonzero masses even in chirally symmetric world due to 
gluons. The same will also perfectly hold true for tech- 
nipions which would always have nonzero masses due to 
technigluons since technicolour QCD should manifest the 
confinement mechanism similar to our one in usual QCD. 
Therefore, technicolour theories look rather doubtful from 
the point of view of our confinement mechanism. 

7.3 Preon models 

Another cardinal approach to the problem of masses is 
connected with the preon models (see, e.g., [3T] and ref- 
erences therein). Under this approch quarks, leptons and 
gauge vector bosons are suggested to be composed of sta- 
ble spin- 1/2 preons, for example, existing in three flavours 
and being combined according to simple rules. The main 
theoretical objection to preon theories is the mass para- 
dox which arises by virtue of the Heisenberg's uncertainty 
principle. Scattering experiments have shown |12] that 
quarks and leptons are point-like up to the scales of or- 
der 10 -3 fm which corresponds to a preon mass of order 
197 GeV (due to the uncertainty principle) if the preon 
is confined to a box of such a size, i.e. its mass will ap- 
proximately be 0.4 x 10 5 times greater than, e.g., that 
of d-quark. Thus, the preon models are faced with a mass 
paradox: how could quarks or electrons be made of smaller 
particles that would have masses of many orders of mag- 
nitude greater than the fundamental fermion masses? The 
paradox might be resolved by the rather dubious postulate 
about a large binding force between preons cancelling their 
mass-energies. Our confinement mechanism points out the 
more physically acceptable way of overcoming these ob- 
stacles. If the interaction among preons is decribed by a 
QCD-like theory based on, e.g., SU(N)-group with N > 2 
then, according to our results 2,3 j, such theories should 
also manifest confinement to generate masses decribed by 
relations similar to (5) and (9). This signifies that preons 
might possess small masses or be just massless and, as a 
result, mass paradox would be removed. 

8 Concluding remarks 

The results of present paper as well as the ones of [HOEIE 
8,9,10 allow one to speak about the fact that the confine- 
ment mechanism elaborated in [TJ[31[3] gives new possibil- 
ities for considering many old problems of hadronic (me- 
son) physics (such as nonperturbative computation of de- 
cay constants, masses and radii of mesons, chiral symme- 
try breaking and so forth) from the first principles of QCD 



immediately appealing to the quark and gluonic degrees 
of freedom. This is possible because the given mechanism 
is based on the unique family of compatible nonpertur- 
bative solutions for the Dirac- Yang-Mills system directly 
derived from QCD-Lagrangian and, as a result, the ap- 
proach is itself nonperturbative, relativistic from the out- 
set, admits self-consistent nonrelativistic limit and may 
be employed for any meson (quarkonium). Under the cir- 
cumstances the words quark and gluonic degrees of free- 
dom make exact sense: gluons come forward in the form 
of bosonic condensate described by parameters aj , bj , Bj 
from the unique exact solution (3) of the Yang-Mills equa- 
tions while quarks are represented by their current masses 
m q . Though nature of the latter is not yet totally under- 
standable (see section 7) but the confinement mechanism 
under discussion indicates a possible way of overcoming 
this puzzle - quarks might be composed from just massless 
preons whose interaction would be described by a gauge 
theory with confinement mechanism similar to that under 
discussion and the current quark masses might be gener- 
ated along the lines discussed in sections 5 and 7. 

The given paper to a great degree summarizes study- 
ing nonet of light pseudoscalar mesons realized in [5,6,7, 
8,10 within the framework of our approach and we can 
ascertain the fact that, on the whole, this nonet can be 
described from the united point of view of our confine- 
ment mechanism. In line with the above, obviously, one 
should now pass on to vector mesons (p, <fr, lo...) and also 
to the light scalar mesons whose nature has been contro- 
versial over 30 years [52]. As is clear from Section 2 (see 
also Appendices A, B), there exists a large number of rel- 
ativistic bound states in the confining SU(3)-gluonic field 
(3) so all the mentioned mesons can probably correspond 
to some of those states and be described by their own 
sets of parameters a,j, bj, Bj of solution (3). More impor- 
tant task is, however, to explore possible ways to extend 
the approach over baryons, in particular, over nucleons. 
In this situation we shall have to deal with a relativistic 
3-body problem as follows from SQM. It is clear, however, 
that confinement mechanism under discussion (which in 
essence decribes the relativistic 2-body problem) should 
also occupy a fitting place in the 3-body constructions. 
We hope to develop the given direction elsewhere. 



Appendix A 

We here represent some results about eigenspinors of the 
Euclidean Dirac operator on two-sphere S 2 employed in 
the main part of the paper. 

When separating variables in the Dirac equation (4) 
there naturally arises the Euclidean Dirac operator 2?o on 
the unit two-dimensional sphere § 2 and we should know 
its eigenvalues with the corresponding eigenspinors. Such 
a problem also arises in the black hole theory while de- 
scribing the so-called twisted spinors on Schwarzschild and 
Reissner-Nordstrom black holes and it was analysed in 
Rcfs. [3"ll2"3"]. so we can use the results obtained therein for 
our aims. Let us adduce the necessary relations. 
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The eigenvalue equation for corresponding spinors <P It should be noted that square of Dq is 
may look as follows 

_ 2 cos$ 1 1 

©o* = Af. (A.l) 2>o = -4^ 2 + ^3-^ + + (A.8) 

As was discussed in Refs. [23], the natural form of T>q while laplacian on the unit sphere is 
(arising within applications) in local coordinates 0, ip on 

the unit sphere § 2 looks as A 1 a ■ Qa . 1 a 2 a2 , . Q q , 1 c,2 

sm0 sin 1? ^ sm i? v 



ict 2 9^ + W7 3 - — „ (dtp — -(T2C3 cos$ 



sin0 V ^ 2 



(A.9) 

so the relation (A.8) is a particular case of the so-called 
Weitzenbock-Lichnerowicz formulas (see Refs. [25]). Then 
1 cot0 from (A.l) it follows T>\(S> — \ 2 <P and, when using the 

aicr 2 d0 + - — -aia 3 d v H — cricr 2 (A. 2) / p, \ 

sm,y 2 ansatz <2> = P(0)e- im v = [p J e ~" n p i,2 = ^1,2(1?), 

with the ordinary Pauli matrices \ 2 / 

the equation 2?q<£ = A turns into 



-flg - cot 00* + , 2 a 4 + . 2 ai P 
v sin u sin w 

so that CTi2?o = — 2?oO"i. 



The equation (A.l) was explored in Refs. [23J. Spec- / \ 



trum of Z?o consists of the numbers A = ±(Z + 1) with ( ^ — 7 J -P , (A. 10) 

multiplicity 2(1 + 1) of each one, where I = 0, 1, 2, .... Let 

us introduce the number m such that — Z < to < Z + 1 and wherefrom all the above results concerning spectrum of 

the corresponding number m' = to — 1/2 so |rn'| < Z + l/2. V can be derived [23]. 

Then the conforming eigenspinors of operator V are When calculating the functions P^, n , (cos 0) directly, 

to our mind, it is the most convenient to use the integral 

>i \ _ ^ _ C ( p t'-i/2 ± P m'i/2 \ - im > v expression [23] 



* = U2 / = ^ TA = 2" 1 P fc , T? 1 , I 

\^ 2 / z y r m'-l/2 + r m'l/2 y 1 , 

(A-3) f* /n ,(coBtf)= 1 Uk-m')\(k + m>)\ 

with the coefficient C = ,/^±± and Zc = Z + 1/2. These 2?r V ( fc ~ n '^ k + n ') ! 

spinors form an orthonormal basis in L 2 ,^ 2 ) and are sub- 27r . x 

ject to the normalization condition / [ e «¥>/ 2 cos |_ i e - i, p/ 2 s j n _ 



2 2 



7T 27T 



f ***dn = J Jm 2 + m 2 )^d^ = i. (A.4) f ie w 2sin !! + e ^/ 2co ^y + "' e w,^ (A . U) 

n n V 2 2 / 





Further, owing to the relation a{D a = -V Q a x we, obvi- and the symmetry relations (z = cos0) 
ously, have 

( T 1 <P TA = $ ±x . (A.5) Pm'M = P n>m-(z)> P m>,-n'( z ) = P -m>,n'( z ), 

As to functions P£ , n , (cos 0) eee P£, | n , (cos <&) then they k = k 

can be chosen by miscellaneous ways, for instance, as fol- m'n'K > -m',-n'\ )■> 

lows (see, e. g, Ref. [23J) .afc-W-2n' ™fc , n im 

^n'H = « p L',-rAz)- (A.12) 

fc , . _ ._,„'_„' \{k — m')\(k — n')\ /l + cos^\ TTi particular 
P m , n ,(costf,-z V (k + miy.(k + n>y. U-costfJ 

V P,\(z) = cos 2 * (0/2), P£_ fc (z) = z 2fc sin 2fe (0/2), 

A (fc + j)!i 2 J / 1-COS0 V ft 

X . Y< V 2 J ^' sin^, 

j— max(m'y) ku \ / 2 fc' 

(A.6) 

with the orthogonality relation at m',n' fixed fe fe , / (2fe)! 



(k - n')l(k + n')\ 



P*V(cos0)P^,„,(cos0)sin0d0= — ^-4fc/ . (A.7) 
mnV ;m " V ^ 2fc + l v ; sin fc -" (z9/2)cos fc+n (0/2). (A.13) 
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Eigenspinors with A = ±1, ±2 



Appendix B 



If A = ±{l + 1) = ±1 then I = and from (A.3) it follows 
that k = 1 + 1/2 = 1/2, \m'\ < 1/2 and we need the 

1/2 

functions P m , _ tl , 2 that are easily evaluated with the help 
of (A.11)-(A.13) so the eigenspinors for A = —1 are 



C/cosf+zsinf ]eiv/2 



-iip/2 



Cf cosf + z , 
2 I -cos^+isin^ ' 



(A.14) 



2 1 " 2 

while for A = 1 the conforming spinors are 



O / cos - - t sin ^ ^ p i V /2 



= cos 2 - t sin 5 
2 I cos£ +z'sin£ 



# = 



L> f — cos ^ + z sin 



2 \ cos f + i sin £ 



2 I 



(A.15) 



with the coefficient C = ^/1/(2tt). 

It is clear that (A.14)-(A.15) can be rewritten in the 
form 

A = -1 : # = ? f e l ^1 e^ 2 



2 V e~ l 2 



or 



C 



* = 2 ( 4-4 1 e ^ /2 



or 



A = 1 : <Z> 



$ = — 



C 



e „ 2 ] e i V /2 



-iip/2 



(A.16) 



so the relation (A. 5) is easily verified at A = ±1. 

In studying vector mesons and excited states of heavy 
quarkonia eigenspinors with A = ±2 may also be useful. 
Then k = I + 1/2 = 3/2, |m'| < 3/2 and we need the 

3 /2 

functions P m , ±1 , 2 that can be evaluated with the help of 
(A.11)-(A.13). Computation gives rise to 



P. 



3/2 



V3 . 



3/2,-1/2 



p; 



.3/2 



sin sin • 



i) 



P 



3/2 

■3/2,1/2' 



3/2,1/2 



• V3 . , 

z — — sin w cos — = P 



3/2 



4^ = 4 (sinf-3sin^l p" 



2 * -3/2,-1/2' 

1/2,1/2' 



p; 



3/2 



1 



1/2,1/2 4 



I? 



1 ( cos " + 3 cos rtf ) = p!{ 2 (A.17) 



/2 -1/2 



We here adduce the explicit form for the radial parts of 
meson wave functions from (6). At rij = they are given 

by 



F jl =C j P j r^e 



-/V 



F 



.72 



/r,c;,r" e 1 ' 1 • 



1 - 



Y, 



(B.l) 



while at rij > they are given by 



F jl =C j P J r»'e 



Y, 



1 



Zi 



F j2 =iC j Q j r^e-^ r x 



2a j +1, 



r J L n,-l 



(B.2) 



with the Laguerre polynomials L^(rj), rj — 2/3j*r, /?,• 

)j Fj 



/i 2 , — to 2 + g 2 b 2 at j = 1,2,3 with & 3 



-(61+&2). P 



g6j + Qj //:: Uj, Y,j = PjQjaj + (P? - Qfjgaj/2, 
Zj = PjQ 3 A 3 + (P? + Qfjga^/2 with a 3 = -(a x + a 2 ), 



.4 



A, 



^ 2 



.gPj with P 3 - V „ A , - y a -j, 

while Xj = ±(lj + 1) are the eigenvalues of Euclidean 
Dirac operator T>q on unit two-sphere with lj — 0, 1, 2, ... 
(see Appendix A) and quantum numbers rij — 0,1,2,... 
are defined by the relations 



gbjZj 



(B.3) 



which entails the spectrum (5). Further, Cj of (B.1)-(B.2) 
should be determined from the normalization condition 



(|F;i| 2 



|F, 2 | 2 



)dr 



(B.4) 



and according to (A.3) this entails eigenspinors with A = 2 
in the form 

2 2 \ e l 2 J 8 \3e'2 +e '2 

C {-3e- i ™-ef\ e - iv /2^C. : /3 s . n{f (- e -^ ^ , ( . ., 



As a consequence, we shall gain that in (4) tf'j s L\{ 
at any t e 1 and, accordingly, & = (^1,^2,^3) may de- 
scribe relativistic bound states in the field (3) with the 
energy spectrum (5). As is clear from (B.3) at rij = we 
have gbj/(3j = Yj/Zj so the radial parts of (B.l) can be 
rewritten as 



Fji = Cjl','-"' 



ft 



3e 1 ^ + e~ l % J ~ '22 \ e'? 

(A.18) 

with C = 1/\/tt, while eigenspinors with A = — 2 are ob- 
tained in accordance with relation (A. 5). 



More details can be found in Refs. [HE] 



9h_ 

ft 



(B.5) 
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Appendix C 

The facts adduced here have been obained in Refs. [21 
[3] and we concisely give them only for completeness of 
discussion in Section 2. 

To specify the question, let us note that in general the 
Yang-Mills equations on a manifold M can be written as 



d * F = g{*F A A — A A *F) , 



(C.l) 



where a gluonic held A = A^dx^ 



AfX n dx^ \X n are the 



known Gell-Mann matrices, /i = t,r,d,p (in the case of 
spherical coordinates), a — 1, ...,8], the curvature matrix 
(held strentgh) F = dA + gA A A = F^ v \ a da^ A dx v with 
exterior differential d and the Cartan's (exterior) product 
A, while * means the Hodge star operator conforming to 
a metric on manifold under consideration, g is a gauge 
coupling constant. 

The most important case of M is Minkowski spacetime 
and we are interested in the confining solutions A of the 
SU(3)-Yang-Mills equations. The confining solutions were 
defined in Ref. [T] as the spherically symmetric solutions 
of the Yang-Mills equations (1) containing only the com- 
ponents of the SU(3)-field which are Coulomb-like or lin- 
ear in r. Additionally we impose the Lorentz condition on 
the sought solutions. The latter condition is necessary for 
quantizing the gauge fields consistently within the frame- 
work of perturbation theory (see, e. g. Ref. [25]), so we 
should impose the given condition that can be written in 
the form div(A) = 0, where the divergence of the Lie al- 
gebra valued 1-form A = A^dx* 1 — A a X a dx^ is defined by 



the relation (see, e. g., Refs. [27] ) 

div(A) = -Ld,(V6g^A„) 



(C.2) 



It should be emphasized that, from the physical point of 
view, the Lorentz condition reflects the fact of transver- 
sality for gluons that arise as quanta of SU(3)-Yang-Mills 
field when quantizing the latter (see, e. g., Ref. [2"S]). 

We shall use the Hodge star operator action on the ba- 
sis differential 2-forms on Minkowski spacetime with local 
coordinates t, r, ip in the form 

*(dt A dr) — — r 2 siwddd A dp , *(dt A d-d) = sin i?dr A dip , 

*(dt A dp) = ^—zdr A d-d , *(dr A dd) = sintM A dp , 

smw 

*(drAdip) = — -?—-dtAd'&,*{d'&Ad<p) = — | — -dtAdr , 



sint? 



r 2 sin i9 



(C.3) 

so that on 2-forms * 2 = —1. More details about the Hodge 
star operator can be found in [37] . 

The most general ansatz for a spherically symmetric 
solution is A = A t (r)dt + A r (r)dr + A^{r)d-d + A v (r)dp. 
But then the Lorentz condition (C.2) for the given ansatz 
gives rise to 

sin ■&d r (r 2 A r ) + ^(sintfAj) = 0, 



which yields A r = Q — J A$(r)dr with a constant 
matrix C . But the confining solutions should be spheri- 
cally symmetric and contain only the components which 
are Coulomb-like or linear in r, so one should put C = 
A$(r) = 0. Consequently, the ansatz A = A t (r)dt+A v (r)d(p 
is the most general spherically symmetric one. 

For the latter ansatz we have F = dA + gA A A = 
—d r A t dt A dr + d r A v dr A dp + g[A t , A v ]dt A dip, where [•, •] 
signifies matrix commutator. 

Then, according to (C.3), we obtain 



*F = (r 2 smti)d r A t dd A dip 



1 



sini? 



d r A v dt A dd- 



sint? 



[A t ,A v ]dr Add , 



(C.4) 



which entails 



d*F = sindd r (r 2 d r A t )drAddAdpA -d 2 A v dtAdr Adtf, 

sin v 

(C.5) 



while 



*FAA-AA*F = r 2 smd\d r A t , A 



10 1 



[SrA^, ) dtAddAdp 



([[A t ,A v ],A t ] dtAdr A d-d + [[A t , A v ], A v ] dr A d-d A dp) . 

(C.6) 



Under the circumstances the Yang-Mills equations (C.l) 
are tantamount to the conditions 



d r (r 2 d r A t ) 
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2 „ Q ll-A-t, A v ], A v ], 
d 2 r A v = -g^At^^At], 



sin^ $ 1 

,2r 



r 2 sm$[d r A t , A t ] 



1 



sini? 



[d r A v ,AJ =0. 



(C.7) 

(C.8) 
(C.9) 



The key equation is (C.7) because the matrices At, A v de- 
pend on merely r and (C.7) can be satisfied only if the 
matrices A t — A^Xa and A v = A^X a belong to the so- 
called Cartan subalgebra of the SU(3)-Lie algebra. Let us 
remind that, by definition, a Cartan subalgebra is a max- 
imal abelian subalgebra in the corresponding Lie algebra, 
i. e., the commutator for any two matrices of the Cartan 
subalgebra is equal to zero (see, e.g., Ref. [35] ). For SU(3)- 
Lie algebra the conforming Cartan subalgebra is generated 
by the Gell-Mann matrices A3, As which are 




1 



10 



,A 8 = -=I0 1 

v d VO -2, 



(CIO) 



Under the situation we should have A t = A\ A 3 + A\ As 
and A v = A 3 V X A + A^X$, then [A t , A v ] = and we obtain 



d r (r 2 d r A t ) =0,d 2 A v = 0, 



(C.ll) 



while (C.9) is identically satisfied and (C.ll) gives rise to 
the solution (3) with real constants Oj , Aj , bj , Bj parametriz- 
ing the solution which proves the uniqueness theorem of 
Section 2 for the SU(3) Yang-Mills equations. 
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The more explicit form of (3) is 

A* = [(a 2 -a 1 )/r + A 1 -A 2 ]/2, 

A 8 t = [Ai +A 2 - (01 + a 2 )/r] V3/2 , 
A% = [(61 - b 2 )r + B x - B 2 }/2, 

Al = \{b 1 + b 2 )r + Bx+B 2 ]^/2 > /2. (C.12) 

Clearly, the obtained results may be extended over all 
SU(A)-groups with N > 2 and even over all semisimple 
compact Lie groups since for them the corresponding Lie 
algebras possess just the only Cartan subalgebra. Also we 
can talk about the compact non-semisimple groups, for 
example, U(iV). In the latter case additionally to Cartan 
subalgebra we have centrum consisting from the matri- 
ces of the form aly (In is the unit matrix N x N) with 
arbitrary constant a. 

The most relevant physical cases are of course U(l)- 
and SU(3)-ones (QED and QCD). In particular, the U(l)- 
case allows us to build the classical model of confinement 
(see Section 2 and Ref. [25]). 

At last, it should also be noted that the nontrivial 
confining solutions obtained exist at any gauge coupling 
constant g, i. e. they are essentially nonperturbative ones. 
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